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Abstract. The nature of non-exponential kinetics in microbial cells inactivation by 
pulsed electric fields (PEF) is discussed. It was demonstrated that possible mechanism 
of non-exponential kinetics can be related to orientational disorder in suspension of 
microbial cells of anisotropic form. A numerical studies of spheroidal cell suspensions 
was carried out. The most pronounced deviations from the exponential kinetics were 
observed for disordered suspensions of prolate spheroids at small electric field strength 
E or at large aspect ratio a. For partially oriented suspensions, efficiency of inactivation 
enhances with increasing of order parameter and field strength. A possibility of the 
PEF-induced orientational ordering in microbial suspensions is discussed. 
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1. Introduction 

Application of high intensity pulsed electric field (PEF) as a nonthermal preservation 
method have been a topic of growing interest over the past decade [fl [2J- PEF 
processing is a promising non-thermal alternative technology and many investigators 
have shown it effectiveness for killing bacteria in liquid aqueous media [HI HI El E]- 
The application of high electric field microsecond pulses (typically 10-50 kV/cm) allows 
inactivation of microorganisms and enzymes [3, [8j El HDl EEI]- However, the mechanism 
of inactivation are not yet fully understood. The important problem is to elucidate how 
kinetics and degree of inactivation depend on the type and geometrical parameters of 
microorganisms, properties of liquid media and their flow regimes, the temperature and 
treatment protocol (electric field strength E, wave forms, pulse duration t iy total time 
of treatment) [T2l[T3]. 

The microbicidal effect of PEF treatment is related to selective damage of the 
biological membrane in microorganism. In an external electric field E, a transmembrane 
potential u m is induced on membrane. When transmembrane potential exceeds some 
threshold value (typically about 0.2-1.0 V), an electric field cause a temporary loss of 
the semipermeability of cell membranes or their electroporation [H1[T5]. The sufficiently 
strong electric field and long time of PEF treatment leads to complete membrane damage 
and cell death [16] . 

The surviving fraction S(t,E), defined as the ratio of the number of undamaged 
microbial cells to the total number of microbial cells, decreases with PEF treatment time 
t. Different empirical models, such as Fermi, Hulsheger, log-log and log-logistic, were 
proposed for description of inactivation kinetics [HI El [HI [19] . Although these models 
are very popular, they have no theoretical justification. The Weibull distribution was 
successfully applied in a number of works for fitting experimental PEF inactivation data 
for a surviving fraction S(t,E) [H EU E21 EH1 EI] 

S(t,E) = eM-(t/r(E)) n W). (1) 

Here, r(E) and n(E) are the time and shape parameters accounting for the effective 
inactivation time and survival curve concavity, respectively [25]. But the Weibull 
model is also empirical and physical meaning of the obtained parameters r(E) and 
n(E) was not elucidated yet. Lebovka and Vorobiev [26] proposed a theoretical model 
for description of the surviving curves of spherically-shaped bacteria with the cell size 
distribution. It was shown that the Weibull model can be successfully applied for fitting 
of the surviving curves during PEF treatment. 

Microbial cells display a variety of shapes and dimensions [271 [28] depending 
on the culture condition and age. Size of cells varies between 0.1 fim and 10 /jm. 
Generally, the following different shapes can occur: near-spherical, ellipsoidal or ovoid 
(cocci), cylindroidal (bacilli), and spiral or comma-like (spirilli). For example, cells of 
Escherichia coli and Salmonella typhi are rod-like and have 0.4-0.6 /xm in diameter and 
2-4 (Ltm in length, cells of Leptospira spp. are very long rod with 0.1 /im in diameter 
and 20 /im in length, cells of Staphylococcus spp are spherical cells with diameter of 
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0.5-1.5 /dm, Sacchoromyces cerevisiae have ellipsoid cells with the principal dimensions 
of 2-8 fim and 3-15 /im, respectively, cells of Klebsiella pneumoniae are ovoid with a 
mean dimension 0.4 /im and Vibrio cholerae have comma-like cells with the principal 
dimensions of 0.5 /im and 1.5-3 /im, respectively [271 [28] . 

In general case, the surviving kinetics may be rather complex. An electroporation 
could be influenced by the aggregation of cells, their arrangement, local cell density, 
local solute concentration, and distribution of local electric field [29], [3Ql EH [32]. The 
killing probability of non-spherical cells depends substantially on spatial orientation and 
changes from cell to cell [9j [33] [31] and can be related with cell diameters, spatial and 
orientational distributions. 

In this work, a theoretical model allowing to describe the survivor curves for 
disordered or partially oriented non-spherical bacteria is formulated. The model predict 
how the cell orientation influences the lifetime of a spheroidal microbial cell exposed by 
PEF. The numerical simulations of surviving kinetics of disordered and partially ordered 
suspensions of microbial cells were done. 

2. Description of the model 

2.1. Transmembrane potential 

In general case, electroporation consists of different stages including the charging of the 
membrane, creation of pores and evolution of pore radii |35j. For a single spherical cell 
under the steady state conditions, the transmembrane potential depends on the angle ip 
between the external field E direction and the radius- vector r on the membrane surface 

u m = 1.5 / RE cos<f. (2) 

Here, R is the cell diameter, and / is a parameter depending on electrophysical 
and dimensional properties of the membrane, cell and surrounding media. In dilute 
suspensions of cells, the value of / is close to 1. 

The value of u m is proportional to the cell radius R. The highest drop of potential 
occurs at the cell poles and decreases to zero at ip — ±7r/2. So, the larger microbial cells 
get killed before smaller ones and the damage probability is maximal at the cell poles. 

If a cell is non-spherical, the transmembrane potential u m becomes more complex 
function of the cell size and geometry, direction of external field and position on the 
membrane surface. The transmembrane potential u m of an arbitrary oriented ellipsoidal 
cell at some point on the membrane surface r(x,y,z) may be calculated from the 
following generalized Schwan equation [38] : 

u m = J2 r^/il-U). (3) 

i=x,y,z 

Here, Lj are the depolarising factors defined by the cell radii R x , R y and R z [39] . 
This approximation works for a membrane with negligibly small conductance and its 
application was discussed extensively in literature [SI EQj HH H21 H3, E] . 
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The depolarizing factor of a prolate spheroid (R z > R x = R y ) along the symmetry 
axis z is (39] 

L. = ^(ln[±f -2e), e = ^a~\ (4) 
and for an oblate spheroid [R z < R x = R y ) it makes 

L z = — (e — Vl — e 2 arcsine), e = Vl — a 2 , (5) 

e 6 

where a = R z /R x is an aspect ratio (major / minor axis). The depolarizing factors in 
x and y directions are defined as 

L x = L y = (1 - L 2 )/2. (6) 

Here, L x = L y = L z = 1/3 for a spherical cell, L x = L y ^ 0.5, L 2 « for a long 
cylinder and L x = L y « 0, L 2 « 1 for a thin disk. 

In general case, the transmembrane potential can be calculated from ([3]), but a 
simpler form of this equation may be considered due to the symmetry of spheroid |34j. 
When the electric field vector lies in a XOZ plane (see figured]), ([3]) can be rewritten as 

u m = xE sin 61/(1 - L x ) + zE cos 0/(1 - £,*), (7) 

where # is an angle between the external field and symmetry axis of spheroid, and x 
and z are coordinates of a point at the spheroid surface. The values of x and z in a 
spheroidal system are defined as [4"5] 



x 



R x \Jl — rf cos 4>,z = R z rj } (8) 

where — l<?y<l,O<0<27r are the spheroidal coordinates. 
Finally, introducing (JHJ) into (GO), we obtain 

_ 1/3 sing cos -K} 1 cos Bar] 
u m = ERa i ( ^— + - - - ), (9) 

where R is a radius of sphere, that has the same volume V as spheroid (V = 4nR 3 /3 = 
A7rR z R 2 j3). 

2.2. Lifetime of a microbial cell exposed by PEF 

The lifetime of a membrane in some point at the spheroid surface can be estimated on 
the basis of the transient aqueous pore model 



r KUj 2 /kT r ) 

T{u m (6, 77, 4>)) = Too exp — — — 10 

1 + {u m {9,r],<p)/u ) 2 

Here, is the parameter (r — > in the limit of very high electric fields), to and 
7 are the line and surface tensions of a membrane, respectively, k = 1.381T0 -23 J/K is 
the Boltzmann constant, T is the absolute temperature, u a = / (C m (e w /e m — 1)) is 
the voltage parameter (the dimension of u Q is Volts), C m is the specific capacitance of a 
membrane, e w ,e m are the relative dielectric permittivities of the aqueous phase and of 
the membrane, respectively. 
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The lifetime of a spheroidal cell r c depends on the angle 9 between electric field 
direction E and the symmetry axis Z of a spheroid. It can be estimated by a averaging 
of T~ l {u m ) on the spheroid surface: 

1 } 2 f \/ a2 (l _ V 2 ) + V 2 

where 

i 

A = 2itR 2 x f ^a 2 (l-T] 2 )+r] 2 dr] (12) 



is the surface area of a spheroid 

For a prolate spheroid, the surface area is 



a arcsin y/1 — a 2 



A = 2kI%(1+ j==z } ' (13) 

and for a oblate spheroid it makes 

^ r.2/ n a arcsin hy/ a~ 2 - 1 . 

A = 2irRl(l + r^—r )• ( 14 ) 

\J a 1 — 1 

2.3. Surviving probability during a PEF treatment 

A surviving probability of a single spheroid with the angle 9 of its principal axis relative 
to the external electric field E is defined as 

S(t,9) = exp(-t/r c (9)). (15) 

Then, a surviving probability S(t) of the whole suspension with spheroids of 
different spatial orientation can be calculated as: 

l 

S(t) = f f(9)exp(-t/r c (9))dcos9, (16) 



where f(9) is an angular orientational distribution function of spheroids. 

For randomly oriented spheroids f{6) = 1/2. For partially oriented spheroids it is 
useful to introduce an order parameter Q defined as [46J: 

1 ) 

Q = -J f(9)(3cos 2 9- l)dcos9. (17) 

-i 

For perfectly oriented suspension, when all spheroids are completely aligned, Q = 1 
and for randomly oriented suspension the order parameter is zero, Q = 0. 

Disordered suspensions of anisotropic cells may be oriented by the external electric 
or magnetic fields [HI SHI SSI ED] , or by the fluid flow [HU [52], [9] . In the external electric 
field E the angular orientation distribution function f{9) can be estimated as [53| H6] 



= exp(£7* cos 2 9)d cos 9 ^ 

/ exp(U* cos 2 9)d cos 9 
-l 
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where U* is a dimensionlees electrostatic energy of spheroid in the external field E, 



Here, (3 is the electrical polarizability anisotropy of a particle that depends on 
electrophysical properties of the particle and the outer solution. 

An order parameter Q can be calculated by substitution of ATS]) into (IT7|) [51] 



Proceeding from (I16I) - (I20I) . surviving kinetics versus order parameter Q can be 
calculated for partially ordered microbial cells. 

2.4- Details of numerical calculations 

For t c (9) evaluation using (l9l)-( |T2|) . the double integration was done using Simpson's 
quadrature rule. The accuracy of numerical integration was better than 10~ 6 . The 
voltage scale parameter in (fTDl was estimated as u Q ~ 0.17V from data obtained by 
Lebedeva [55] for the general lipid membranes (u ~ 1.69 * 10~ n N, 7^2* 10~ 3 N/m, 
e w Rj 80, e m ps 2, C m ~ 3.5* 10~ 3 F/m 2 at T = 298K). The time scale parameter was put 
as ps 3.7* 10~ 7 s [55]. It is useful to use in calculations a dimensionless reduced field 
intensity defined as E* = E/E Q , where E Q = 2u Q /3R, and R is an equivolume radius 
of a spheroid. Note that at 2R ps l//m, and w ps 0.17V, E a ps* 2.27kV/cm. All these 
parameters were used calculation for estimation purposes. 

The calculated dependencies of t c (9) were used for numerical calculation of the 
surviving kinetics from (TIB"]) at different values of order parameters Q ( (120]) ). 

3. Results and discussion 

3.1. Lifetime of a spheroidal microbial cell 

Figure [2] and figure [3] present some examples of the calculated relative lifetime r c / 
versus reduced field intensity E* for a prolate (|2J) and oblate spheroids at different 
values of angle 9. A prolate spheroid in external electric field was more stable at 9 = 90° 
and less stable at 9 = 0° than a spherical cell of the same volume, but an oblate spheroid 
was always less stable electrically than a spherical cell of the same volume. 

The relative lifetime Tc/t^ versus angle 9 for different aspect ratio a for a prolate 
and oblate spheroids at E*=10 are presented in figure H] and figure [5j 

For a prolate spheroid, the value of Tc/tqc considerably increases with angle 9 
increase and it was a minimum for a cell aligned along the applied field E (figure H]). This 
result is in accordance with experimental observations of [34] . who reported minimum 
electropermeabilization for the cells aligned along to the electric field direction. 

For an oblate spheroid, the value of Tc/t^ was smaller than for a spherical cell of 
the same volume and the value of t c Itqq decreases with angle 9 increase (figure [5]). It 



U* = l3E 2 /(2kT). 



(19) 




(20) 
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was maximal for = 0, but Tc/t^ dependence versus angle 9 was not so distinct as for 
a prolate spheroid. 

There exist some threshold angle 9 = 9t, at which the curve t c /too{E*) for a 
prolate spheroid is very close to that for a spherical cell of same volume. The cell 
permeabilization was suppressed at 9 > 9 t and at was enhanced 9 < 9 t as compared 
with a spherical cell of the same volume. The higher was the aspect ratio a, the larger 
was the threshold angle 9 t . For an oblate spheroid, the threshold angle 9 = 9 t was 
observed only at high aspect ratio a > 0.3. 

3.2. Survivor kinetics of a disordered suspension of cells 

Due to the Brownian motion a random orientational distribution for suspensions of 
microbial cells is typical when field-induced ordering effects are absent. Figure [6] and 
figure [7] show the calculated survivor curves Sit) of disordered suspensions (f(9) = 1/2) 
in (TT6l)). The ideal first order kinetics law ( lT5i) was observed for suspensions of identical 
spherical cells (dashed lines in figure [6] and figure [TJ. The deviations from first order 
kinetics for prolate spheroidal cells became more pronounced with decreasing of the 
electric field intensity E* (figure [6]) or increasing of the aspect ratio a (figure [7j). The 
Tel Too versus 9 dependence was not so pronounced for oblate cells as for prolate cells, 
and no noticeable deviations from the first order kinetics were observed. 

The kinetics Sit) demonstrates that surviving probability in disordered suspension 
was higher for prolate cells and was lower for oblate cells as compared with the surviving 
probability for spherical cells of equivalent volume (figure [7]) . 

The calculated survivor curves Sit) for the prolate cells may be fitted with empirical 
Weibull function ([I]). This model always gives only upward concavity, i.e. n < 1, for 
orientationally disordered suspensions of prolate cells. But the numerically estimated 
shape n and relative time r/r^ parameters were rather sensitive to the upper cutting 
boundary tm^jT^. This fact reflects existence of an intrinsic inconsistency between an 
unknown survival function and Weibull function [26J. 

3. 3. Survivor kinetics of a partially oriented suspension of cells 

Figure [8] and figure [9] show the calculated survivor curves S(t) for partially ordered 
suspensions of spheroidal cells. The surviving kinetics of more disordered suspensions 
(Q — > 0) in the limit of large time it/r^ ^> 1) was obviously controlled by 
the inactivation of cells oriented perpendicular to the applied field, and Sit) ~ 
exp(-t/r(90°)). 

Increase of the order parameter Q results in two different regimes of surviving 
kinetics that correspond to inactivation of cells oriented along the field (fast regime at 
small time Sit) ~ exp(t/r(0°)) and perpendicular to the applied field (slow regime at 
large time S(t) ~ exp(t/r(90°)). The partial contribution of the first (fast) regime to 
inactivation kinetics increases with growth of the order parameter Q (figure [8]). 
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For partially orientationally ordered suspensions with given Q, increase of the 
electric field strength E* causes enhancement of inactivation kinetics and two regimes 
of surviving kinetics are also observed ( figure [H])- The rate of inactivation in the regime 
of slow inactivation is not constant in the limit of large time (t/r^ ^> 1) and increases 
with E* increase. 

The orientational ordering can be induced during the PEF treatment. Because of 
the quadratic dependence (see equations f|T8l) . f|T9|) ) on the electric field strength E, the 
orientational ordering in high pulsed electric fields may be noticeable. The possibility 
of ordering for ellipsoidal or cylindrically shaped microorganisms under the effect of 
external electric fields is discussed in [56J. A rod-shaped tobacco mosaic virus (TMV, 
about 0.018 /im in diameter and 0.3 /im in length [57] ) demonstrates a strong orientation, 
near to complete saturation of the optical birefringence, in the electric field as high as 
~4 kV/cm[5l]. Electrooptical studies of rod-shaped E. coli suspensions [581 E2] shows 
the existance of strong orientational ordering at electric fields of E < 1 kV/cm. 

The orientation electric field induced effects can be roughly estimated using 
(fT9]) . (1201) . Taking the experimental value of the electrical polarizability anisotropy of E. 
coli cells (3 = 410- 27 Fm 2 p], we obtain from (pD that U* = pE 2 /87rkT « 510 3 at 
E = 1 kV/cm and T = 298 K. The corresponding order parameter following from ( 1201) 
is Q ~ 1, e.i., degree of ordering is high. 

But the degree of orientation can depend also on the pulse duration. As it is shown 
schematically in figure [10], the order parameter Q increases with time constant r Q after 
the external electric field is switched on. The order parameter Q decreases to zero with 
another time constant r x when the external field is switched off. The relaxation time r x 
is determined by the Brownian rotation diffusion of the spheroid rotation about x axis 
in absence of electric field. 

According to [60], the rotational diffusion times of a spheroid with respect to the 
symmetry axes x (or y), and z are: 

2 (« 4 - 1 ) 

Tx ' y ~ 3a((2a 2 -l)F-a) TiJ ' ( } 

^r, (22) 



3a(a — F) 

where tr =7ri]R 3 /kT is the rotational diffusion time for a sphere of radius R with the 
same volume spheroid and 



F = F p = ln(a + Va 2 - 1) /Va 2 - 1, (23) 

F = F = arctan(Va- 2 - 1) /VI - a 2 , (24) 

for a prolate and oblate spheroid, respectively. 

Insert in figure [TUl shows a rotation diffusion time r x (along the short axis x) of the 
prolate spheroid versus an aspect ratio a as calculated from ( 12T1) . ( l23l) . For example, 
the rotational diffusion time is of the order of r x ~ls for a E. coli cell with equivolume 
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radius R m 0.64/im and aspect ratio a ~ 2 [55], but it can increase substantially with 
increase of the equivolume radius R or aspect ratio a. 

It is more difficulty to calculate the time r a characterizing the process of ordering 
in the external electric field . The estimations shows [l6]that r Q m t x in a low field, 
when, Q < 0.1. At a very high electric field, when Q — > 1.0, t q ~ t x /(2C/*/15), where 
U* = f3E 2 /(8rrkT) (equation ([19])). 

For rather small E. coli cells from the above estimations, we obtain r « 
t x /(2U*/15) « 10~ 3 s at £ G = 1 kV/cm and r « KT 5 s at £ = 10 kV/cm (T = 298 K). 
So, small cells like those of E. coli may be effectively oriented during the pulse duration 
at PEF treatment with field strength E — 10 kV/cm and pulse duration £j = 10 _5 s. 
This estimation is in accordance with experimental observations (62] showing that the 
24 ms pulse causes the observable orientation of the E. coli cells parallel to the field 
direction only at fields E exceeding 1.25 kV/cm. 

However, the PEF-induced orientation effects can be supressed for larger cells 
or in presence of cells aggregation. The factor of bacterial aggregation is essential, 
because particles of the large colloidal aggregates are bounded and can not be freely 
reoriented during the PEF treatment. The aggregation and colony formation is a typical 
phenomenon in biocolloidal suspensions of bacterial particles, and some kind of bacteria 
(e.g. Bacillus subtilis) exhibit various aggregation patterns [63j El]. Moreover, the 
PEF-induced cell-cell aggregation are also inportant [65J. 

For larger cells, for example, with equivolume radius R of « 5.0/im the rotation 
diffusion time is r x ~ 10 2 — 10 3 s (insert in figure fTUi) . In this case r G t x /(2U* /15) ~ 
10~ 2 — 10~ 3 s at E = 10 kV/cm and ij <C r Q . So, PEF-induced orientational effects 
can be rather small and not contribute into inactivation kinetics for large cells or cell 
aggregates. 

4. Conclusion and outlook 

This study has demonstrated the possibility of non-exponential survaving kinetics of 
microbial inactivation under the PEF treatment, which is believed to be related to the 
orientational disorder existing in a suspension of spheroidal microbial cells. Another 
factor can be related with sizes distribution of microbial cells [26J. Deviations from 
the ideal first order kinetics law (( fl~5i) ) are more pronounced in completely disordered 
suspensions of prolate spheroids at small electric field strength or at large aspect ratio a. 
Efficiency of inactivation enhances with order parameter and field strength increase in 
partially oriented suspensions. In general case the inactivation kinetics can be influenced 
by the concentration of cells, their aggregation and arrangement [30l l3T| [32]. The 
relevance and importance of such effects for explanation of the survival curves observed 
in PEF-inactivation experiments should be studied in future. The accurate description of 
inactivation kinetics requires accounting for the dynamics of bacterial cell reorientation 
in a high electric field during its inactivation. It seems to be important also to find 
correlations between factors that controls aggregations of bacterial cells, PEF protocols 
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and parameters of bacterial inactivation kinetics. 
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Figure legends 




Figure 1. Schematic representation of the problem under consideration. X, Y, Z are 
the local Cartesian coordinates related to the oblate spheroidal microbial cells in the 
external field E (it is supposed here that it is lies in the XOZ plane). Here, 9 is an angle 
between the electric field direction E and symmetry axis of spheroid Z , r(x, y, z) is the 
radius of a membrane surface point, where the transmembrane potential is calculated, 
Rx = R y ,Rz are the cell radii. 
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Figure 2. Relative lifetime Tc/t^ of a prolate (aspect ratio a = 10) cells versus 
reduced field intensity E* at different angles 9 between electric field direction and 
symmetry axis of spheroids. The solid lines show data for a spherical cell with a radius 
equivalent to that of an equivolume spheroid. 
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Figure 3. Relative lifetime t c /too of a oblate (aspect ratio a = 0.10) cells versus 
reduced field intensity E* at different angles 8 between electric field direction and 
symmetry axis of spheroids. The solid lines show data for a spherical cell with a radius 
equivalent to that of an equivolume spheroid. 
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Figure 4. Relative lifetime Tc/tqc of a prolate spheroidal cell versus angle 9 for 
different aspect ratio a. The calculation were done at the given value of reduced 
electric field intensity E* = 10 that corresponds to E w 22.7 kV/cm for an equivolume 
radius of spheroid 2R « 1/im, and voltage parameter u Q ss 0.17V [55 . Arrows show 
threshold angles 0f 




Figure 5. Relative lifetime t c /toc of a oblate spheroidal cell versus angle 9 for different 
aspect ratio a. The calculation were done at the same condition as for data in figure |U 
Arrows show threshold angles 8 t . 




Figure 6. Survivor curves S(t/Too) for orientationally disordered spheroids at given 
a = 10 and different E* . The dashed lines show data for spherical cell cquivolumc with 
spheroids. 




Figure 7. Survivor curves S , (i/r 00 ) for orientationally disordered spheroids at given 
E*~10 and different a . The dashed lines show data for spherical cell equivolume with 
spheroids. 




Figure 8. Survivor curves 5(t/roo) for partially ordered suspensions of prolated 
spheroids at a = 10, E* = 10 and different Q values. The dashed lines show the 
slopes that correspond to orientation with 8=0 and 9 =90°. 




Figure 9. Survivor curves 5(t/roo) for partially ordered suspensions of prolated 
spheroids at a = 10, Q = 0.65 and different E* values . The dashed lines show 
the slopes that correspond to orientation with 9=0 and 9 =90°. 
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Figure 10. Changes of order parameter Q versus time t in the external electric 
field. Here, t x is the diffusion time for spheroid rotation about x axis, r Q is the same 
diffusion time in the presence of the external field. Insert shows t x versus aspect ratio 
a for prolate spheroids, estimated from $2$ and flSSJ), T — 298 K, i] — 8.911CT 4 Pa s 
(water viscosity [61]), R is the radius of a sphere with the same volume as spheroid. 



